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ABSTRACT 

Context. The standard unified scheme of active galactic nuclei requires the presence of high column densities of gas and 
dust potentially obscuring the central engine. So far, few direct subarcsecond resolution studies of this material have 
been performed toward radio galaxies. 

Aims. The goal of this paper is to elucidate the nuclear environment of the prototypical X-shaped Fanaroff- Riley type 
II radio galaxy 3C403, the only powerful radio galaxy known to host an H2O megamaser. 

Methods. Very Large Array A-array and single-dish Green Bank and Effelsberg 1.3 cm measurements were performed 
to locate and monitor the water maser emission. Very Long Baseline Interferometry 6 cm continuum observations were 
taken to analyze the spatial structure of the nuclear environment at even smaller scales, while the CO J=l— and 2-1 
transitions were observed with the IRAM 30-m telescope to search for thermal emission from a spatially extended, 
moderately dense gas component. 

Results. Positions of the H2O maser features and the continuum emission from the core coincide within 5 mas (5.5 pc). 
Intensities of the two main maser components with (isotropic) luminosities sometimes surpassing 1000 L0 appear to be 
anti-correlated, with typical timescales for strong variations of one year. If the variations are intrinsic to the cloud(s), 
the implied angular source size would be JS0.3 mas and the brightness temperature <^5xl0 8 K. The VLBI continuum 
observations support a scenario where a nuclear core, represented by the dominant central radio continuum component, 
is accompanied by a jet and counterjet, directed toward the western and eastern large scale lobes of the galaxy. CO 
remains undetected, providing a maximum scale size of ~50pc x (500 K/Tb) 1 ^ 2 , with Tb denoting the brightness 
temperature of the CO J— 1-0 line. Possible scenarios that could produce the observed maser emission are outlined. 
Adopting a mass of several 10 8 for the nuclear engine, the observed maser features can only be interpreted in terms of an 
accretion disk as in NGC4258, if they solely represent the systemic velocity components. The receding and approaching 
parts of the putative maser disk are, however, not seen and a secular velocity drift of the observed features is not 
(yet) apparent. Most likely, the two main maser components mark shocked molecular gas interacting with the nuclear 
jets. The X-shaped morphology of the radio galaxy may point at a binary nuclear engine. This possibility, greatly 
complicating the nuclear environment of 3C403, should motivate a number of worthwhile follow-up studies. 

Key words. Galaxies: individual: 3C 403 - Galaxies: active - Galaxies: ISM - masers - Radio lines: ISM - Radio lines: 
galaxies 



1. Introduction 

The physical conditions in active galactic nuclei (AGN) are 
unique in the cosmos. Stellar and gas densities are excep- 
tionally large and enormous amounts of energy and angu- 
lar momentum are released as the material accretes onto 
the supermassive nuclear engine. Performing studies of the 
structure, kinematics and excitation of this material is the 
sole mean available to investigate (super)massive ultracom- 
pact objects. 

Because of obscuration and crowding, only a few molec- 
ular tracers can provide direct information on the in- 
nermost regions of active galaxies. The H2O A=1.3cm 
maser line, profiting from a long wavelength (no obscura- 
tion by dust grains), from its enormous luminosity (up to 
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10 4 Lq (isotropically), corresponding to ^10 53 photons/s; 
see Koekemoer et al. 119951 Barvainis & Antonucci 2005) 
and the availability of Very Long Baseline Interferometry 
(VLBI) networks, is a particularly suitable tracer of this en- 
vironment. The line samples gas with kinetic temperatures 
of several 100K and densities of ^10 8 cm~ 3 (e.g., Kylafis 
& Norman [LWlfTMTD . 

So-called " disk-masers" as in NGC 4258 allow us to map 
nuclear accretion disks (e.g. Greenhill et al. 119951 Miyoshi 
et al. 119951 Herrnstein et al. l2005)l . Analysis of the dynamics 
of observed Keplerian disks leads to estimates of the mass of 
the nuclear engine and to a determination of the distance to 
the galaxy, independent of any standard candles. So-called 
"jet-masers" as in Mrk 348 (Peck et al. I2U03[) are associated 
with the pc-scale jets and provide estimates, through rever- 
beration mapping, of the speed of the material in the jet. A 
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further class of masers is associated with nuclear outflows 
(Circinus; Greenhill et al. [2"0"0"3"|) . 

Measuring the proper motion of galaxies has been an- 
other long-standing problem. The solution is to find suf- 
ficiently luminous ultracompact sources that can be ob- 
served with VLBI techniques involving phase referenc- 
ing. Extragalactic H2O masers are suitable targets and 
could provide, through their measured radial velocities and 
proper motions, three dimensional velocity vectors in space 
(Brunthaler et al. 120051 I2007j) . not only for the masers 
themselves but also for the nuclei of their parent galaxies. 

Based on the assumption that the water masers are lo- 
cated in circumnuclear tori, (often) aligned with the nucleus 
and amplifying its continuum, linearly for saturated masers 
and exponentially for unsaturated ones, one should expect 
a high detection rate in radio-loud galaxies (in particular 
those with the radio axis close to the plane of the sky). 
However, water megamasers have been found so far only in 
low-radio power AGN (Pighz < 10 22 " 23 WHz" 1 ), mostly 
in galaxies classified as Sy 2s or LINERs (Braatz et al. 1996, 
[15971) . 

A first systematic search for the A=1.3cm line in ra- 
dio galaxies was performed by Henkel et al. 0998). No 
maser was detected in a sample of ~50 Fanaroff-Riley type 
I (FR I) galaxies. Subsequently, targets belonging to the 
more powerful Fanaroff-Riley type II (FR II) class were 
observed by two groups (Tarchi et al., Lara et al., unpub- 
lished) again with negative results. 

According to the current paradigm for radio-loud AGN, 
FR II galaxies host type 2 AGN with narrow optical emis- 
sion lines and prominent, extended nuclear jets. Since the 
jets are ejected at a small angle with respect to the plane 
of the sky, FR II galaxies should contain a nuclear torus 
being viewed approximately edge-on (for a review, see e.g. 
Urry & Padovani ll995[) . So far, however, the presence and 
nature of these hypothesized tori has remained an enigma. 
The nuclear region of the prototypical FR II galaxy Cyg A 
shows evidence for dust and obscuration with gas column 
densities up to several 10 23 cm -2 . Morphology and physical 
parameters of the nuclear region as well as the mechanism 
triggering nuclear activity are, however, not yet fully un- 
derstood (e.g., Conway & Blanco [T9951 Carilli & Barthel 
Fuente et al. 200f]} Bellamy & Tadhunter [HE]) . 

Very recent modifications to the standard unified 
scheme (hereafter SUS) state that FR I galaxies (or at least 
the majority of them) seem to lack a geometrically and 
optically thick molecular torus. This scenario is now sup- 
ported by an increasing number of studies (e.g., Chiaberge 
et al. lTM9l Perlman et al.[2"00"Tl Verdoes Kleijn et al. 1201121 
Whysong & Antonucci 2004) . In order to account for the ab- 
sence of broad emission lines in narrow-lined FR lis (within 
the framework of the SUS believed to represent the parent 
population of radio-loud quasars) a geometrically and opti- 
cally thick obscuring layer is required instead (Barthel ll989l 
and references therein). 

Recently, a luminous (Lh 2 o^10 3 L q ) H 2 mcgamaser 
has been detected in 3C403 (Tarchi et al. 2003; hereafter 
THC), an FR II galaxy at roughly the same distance as 
Cyg A. 3C403 is the first and so far only powerful radio 
galaxy with a detected H 2 maser. 

The kpc radio morphology of 3C 403 is typical for X- 
shaped radio galaxies. The cause of the "X" radio shape 
is still debated. Arguments exist that support the origin 
of this peculiar morphology as due to a sudden change in 



the jet direction. Merritt & Eckers (j2002[) have argued that 
the sudden change in the jet direction reflects a sudden (< 
few yr) change in the black hole spin axis due to a black 
hole merger (see also Sect. 5.2.2). Alternative models imply 
instead slow precession of the jet axis or lobe backflow (see 
Dennett-Thorpe et al. [20021 Schoenmakers et al. 12000]) . 

The Very Large Array (VLA) A=3.6cm maps of Black 
et al. (|T992l hereaf ter BBL; their Fig. 13) and Dennett- 
Thorpe et al. (|2"0"0"2l hereafter DSL; their Fig. 1) show two 
bright radio lobes towards the east and west with hot-spots 
and two weaker, longer and broader wings to the north- 
west and south-east. The two bright lobes are caused by 
jets containing a number of prominent knots. Adopting a 
W MAP cosmology (H =71, fi M =0.27, ft A =0.73; Spergel et 
al. I2003p . the redshift of z=0.059 leads to an angular size 
distance of Da ~ 230 Mpc, at which 1 mas corresponds to 
~l.lpc. 

In this paper, we present results from a follow-up VLA 
A-array observation of the maser emission and the outcome 
of a search for CO emission using the 30-m IRAM telescope 
at Pico Veleta. We also discuss monitoring observations of 
the maser line spanning two years. Furthermore, to obtain a 
more realistic picture of the innermost regions of 3C 403 and 
to find out whether an accretion disk scenario is plausible, 
we mapped the 6 cm radio continuum emission with parsec 
scale resolution using the European VLBI Network (EVN) . 

2. Observations and data reduction 

Effelsberg 3C 403 was observed in the 616 — 5 2 3 transi- 
tion of H 2 (rest frequency: 22.23508 GHz) with the 100-m 
telescope of the MPIfR at Effelsber£| at four epochs be- 
tween January 2003 and April 2004. The beam width was 
40". The observations were made with a 18-26 GHz HEMT 
receiver that was sensitive to the two orthogonal linear 
polarizations. A dual beam switching mode with a beam 
throw of 2' and a switching frequency of ^1 Hz was chosen. 
The Effelsberg AK90 spectrometer was used as backend. 
System temperatures, including atmospheric contributions, 
were ~40-60K on an antenna temperature scale (T£). The 
beam efficiency was 77b ~0. 5. Flux calibration was obtained 
by measuring W3 (OH) (see Mauersberger et al. H988|) . Gain 
variations of the telescope as a function of elevation were 
taken into account (see e.g., Eq. 1 of Gallimore et al. ()2001|) . 
The pointing accuracy was better than 10". All data were 
reduced using standard procedures of the GILDAS software 
package ( [http://www.iram.fr/IRAMFR/GS/gildas.html ). 

GBT Observations were made with the Green Bank 
Telescope (GBT) of the NRAC0 during six sessions be- 
tween 2003 October 20, and 2005 February 15. We used 
the 18-22 GHz K-band receiver, which has two beams at 
a fixed separation of 3' in azimuth. The GBT beamwidth 
is 35" at the red shifted maser frequency of 21 GHz, and 
pointing uncertainties were better than 10". The data were 
taken in total power mode, and the telescope was nodded 
between two positions on the sky such that the source was 
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always in one of the two beams during integration. We used 
a nod cycle of 2 minutes per position. 

The spectrometer was configured with two 200 MHz 
bandpasses, one centered on the systemic velocity of 3C 403 
and the second redshifted by 180 MHz. No emission was de- 
tected beyond what is shown in Fig.[TJ The zenith system 
temperature was between 35 and 45 K. Atmospheric opacity 
was estimated using system temperature and weather data, 
and ranged from 0.04 to 0.06 at the zenith. We reduced the 
data using GBTIDL. 

VLA Water maser emission was observed at the red shifted 
maser frequency of 21 GHz on July 23 and 25, 2003, with 
the Very Large Array (VLA) of the NRAO 2 in its A con- 
figuration for an 8 hour total on-source time. We observed 
with two 12.5 MHz IFs centered on the two maser features 
detected by Tarchi et al. (|2003l see the upper spectrum of 
Fig.[T]). Each IF was split into 32 channels providing a chan- 
nel spacing of ~ 6 kms -1 . The source 1331+305 (2.66 Jy) 
was used as flux and bandpass calibrator. The point source 
1950+081, at a distance of 5?6 from 3C403 (there was 
no suitable closer source), was used for phase calibration. 
The two IFs were calibrated separately and then combined 
(AIPS task: UVGLU) yielding a total velocity coverage of 
-360 kms" 1 . The data were Fourier-transformed using nat- 
ural weighting to create a 256 x 256 x 64 data cube. The ra- 
dio continuum was subtracted using the AIPS task UVLSF. 
This task fits a straight line to the real and imaginary parts 
of selected channels and subtracts the fitted baseline from 
the spectrum, optionally flagging data having excess resid- 
uals. In addition, it provides the continuum as a UV data 
set, which has been used to create continuum maps of the 
galaxy. The restoring beam is O'.'l x O'.'l. The rms noise 
per channel becomes 0.6 mJy, consistent with the expected 
thermal noise. 

EVN 3C403 was observed on May 20, 2004, with the 
European VLBI Network (EVN) at 6 cm. Eight 8 MHz 
bands, each at right and left circular polarization, were em- 
ployed. The total observing time was 10.5 hours. We used 
4C +02.49 as a phase-reference source and switched every 
2 minutes between the two sources. The initial calibration 
was performed with the AIPS package. A priori amplitude 
calibration was applied using system temperature measure- 
ments and standard gain curves. Fringes were found in the 
3C 403 data itself on all baselines, so we used this source as 
phase-reference and applied the solutions also to 4C +02.49. 
The data were self-calibrated and mapped using the soft- 
ware package DIFMAP (Shepherd et al. H5M|) . We started 
with phase-only self-calibration and later included phase- 
amplitude self-calibration with solution intervals slowly de- 
creasing down to one minute. 

Pico Veleta 3C 403 was observed in a search for CO emis- 
sion on November 6, 2004, using the IRAIvfl 30-m telescope 
at Pico Veleta, Spain. Using the A/B receiver combination 
we searched simultaneously for CO (1-0) and CO (2-1) with 
the A/B 100 receivers tuned to 108.85 GHz and the A/B 230 
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4 IRAM is supported by INSU/CNRS (France), the MPG 
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receivers tuned to 217.69 GHz. Beamwidths were 23" and 
13", respectively. Spectra were obtained using the wobbler 
switch technique with a beam throw of 60" and a switch- 
ing frequency of 0.5 Hz. The data were recorded using the 
1 MHz (512 channels, 1MHz channel spacing) and 4 MHz 
(256 channels, 4 MHz channel spacing) filterbanks for the 
109 and 218 GHz observations, respectively. The effective 
velocity coverage is 1400 km/s at both frequencies. Typical 
system temperatures were 160 and 500 K (Tt). The fast 
switching procedure and the good quality of the resulting 
baselines also provide a rough measure of continuum levels 
(see Sect. 4.1.3). 3C403 was observed for about 70 minutes 
(on+off) leading to rms noise levels of 4.2 and 12.5 mJy 
after smoothing the 109 and 218 GHz spectral bands to 
channel widths of 44 and 50 kms -1 . 

3. The systemic velocity of 3C 403 

For the following sections it is of fundamental relevance to 
assess the accuracy of the reported systemic velocity V sys 
of 3C403. Optical emission line measurements indicate V sys 
= 17688 kms" 1 (z = 0.059) with no uncertainty value pro- 
vided (Spinrad et al. 119851 and references therein). Previous 
studies indicate that recessional velocities of galaxies de- 
rived from optical emission lines tend to have a systematic 
uncertainty, often being smaller than systemic velocities de- 
rived from Hi (e.g. Mirabel & Wilson [19841 Morganti et al. 
I200ip . This difference is attributed to a net outflow of gas in 
the optically unobscured front side of narrow emission line 
regions (e.g. Mirabel & Wilson 11984)) . Unfortunately, no 
neutral hydrogen has been so far detected in 3C 403 (see 
Sect. 5.1.2) to perform a comparison with the optically- 
derived recessional velocity. Baum et al. (|1990p . while not 
providing a recessional velocity, derived on the basis of [Oi] , 
[Nil] , Ha and [Sn] data a rotation curve of 3C 403 that in- 
dicates no strong anomalies. From this we conclude that 
the uncertainty in V sys , as derived by Spinrad et al. 1985, 
should not exceed lOQkms -1 . 

4. Results 

4.1, Spectral line emission 
4.1.1. Maser monitoring 

Figure[T] shows the water maser spectra from 3C403 at 
eleven epochs from January 2003 to February 2005, ob- 
served with the Effelsberg telescope, the GBT, and the VLA 
A-array. Past Effelsberg unpublished data from December 
1997 showed no sign of emission at a 3cr level of 60 mJy. 
Fig. [5] displays our most sensitive spectrum, that obtained 
with the VLA indicating the possible presence of three fea- 
tures that, for the sake of convenience, we indicate as com- 
ponent I, II, and III. 

The most obvious change that can be deduced from 
Fig. [T] is that, in the early spectra, emission is mostly 
present on the red-shifted side (component I), while later 
most of the emission (component II) is blue-shifted w.r.t. 
the systemic velocity. In January and March 2003, com- 
ponent I had a peak line strength of 20-25 mJy, becom- 
ing weaker at subsequent epochs and remaining undetected 
in the most recent spectra (<5mJy). Component II corre- 
spondingly increases its flux from —4 to 20 mJy before sig- 
nificantly fading in the last two spectra. From March 2003 
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Fig. 1. Mascr lines in 3C 403 observed with Effelsberg (four 
epochs), the GBT (6 epochs), and the VLA (one epoch). 
The spectrum from the first epoch (upper panel) was al- 
ready published in THC. The channel spacing is 312 kHz 
(= 4.5kms~ 1 ), 195kHz (= 2.8kms- 1 ), and 391kHz (= 
S.Gkms^ 1 ) for the spectra taken with Effelsberg, the GBT, 
and the VLA, respectively. Zero velocity corresponds to a 
frequency of 20996.28 MHz (cz = 17688 km s _1 w.r.t. the 
Local Standard of Rest (LSR)), the galaxy's recessional ve- 
locity according to the NASA/IPAC database (NED) 



to January 2004, component II seems to experience a first 
shift from about -45 to -95kms _1 and a second one back 
to about -60kms _1 in February 2005. Of course, the pos- 
sibility that there is actually no shift and we are, instead, 
witnessing the brightening and fading of different compo- 
nents is also a possibility. Component III, detected at a low 
significant level only in the first three spectra, fades below 
the 5 mjy noise level together with component I. 

While no systematic regular velocity drift involving all 
the feature can be confidently observed, it is remarkable 
that within about 15 months the maser lines swapped their 
intensities. A time scale of one year and peak flux den- 
sities of order 20mJy indicate an angular source size of 




Velocity (km/s) 



Fig. 2. The main water maser features in 3C403 observed 
with the VLA A-array on July 23, 2003 (see Sect. 11X21 and 
Fig. [J). The dashed vertical line marks the nominal systemic 
velocity of the galaxy, V sys = cz = 1 7688 km s -1 . 
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Fig. 3. Close-up of the region around the nucleus of 3C 403. 
Shown as a filled triangle, a cross, and a filled square are 
the variance-weighted average locations with error bars of 
maser components I, II, and III, respectively. Position off- 
sets in (RAj2oooi Decj2ooo) are relative to the 22-GHz con- 
tinuum peak position. 



^0.3 mas and a peak brightness temperature of <;5xl0 8 K, 
if variations are intrinsic to the cloud(s) and are not trig- 
gered by fluctuations of the continuum background. Each 
individual component reaches peak isotropic luminosities 
well in excess of 1000 L Q . So far, only four other maser 
sources were reported to have a similar power, TXS22226- 
184 (Koekemoer et al. fl995| . Mrk 03 4 (He nkel et al. [2005|) . 
J0804+3607 (Barvainis & AntonucciEDQl]), and UGC 5101 
(Zhang et al. I2006j) . A variation by more than a factor of 
two has not yet been reported in such a luminous maser 
component. Less luminous masers with linewidths in ex- 
cess of lOkms" 1 are commonly more stable. 

4.1.2. The VLA view of the H 2 maser 

In all VLA velocity channels (see Fig. [2 for a preliminary 
report, see also Tarchi et al. 200")) the emission peaks at 
the position RA J20 oo = 19 h 52 m 15?80, Dec J2000 = 02° 30' 



24. 2. The nominally estimated error for absolute positions 
for a VLA A-array map is O'.'l when a strong nearby cal- 
ibrator can be used for phase referencing (as done in our 
observations, using the source 1950+081). Since we have 
also used the "fast switching" technique, this error can be 
taken as a safe upper limit. 
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Fig. 4. 30-m IRAM telescope CO(l-O) (solid lines) and 
CO(2-l) (dotted lines) spectra of 3C403. The data were 
smoothed to channel widths of 44 and 50 kms" 1 , respec- 
tively. The zero velocity corresponds to the recessional ve- 
locity of the galaxy (V sys — 17688 kms -1 ). 



The total integrated intensity (moment-0; in 
mJy kms -1 beam" 1 ) VLA plots of components I— III 
show that they are unresolved at the O'.'l resolution and 
coincident in position with the radio continuum nucleus of 
the galaxy (see also Sect. l4.2.Tj) . The alignment between line 
and continuum peak positions derived from the same data 
is only limited by the signal-to-noise ratios (SNRs) and is 
estimated to be cr reJ = ^(<V 2 • SNR t ) 2 + (0 c /2 • SNR C ) 2 
= ^(0.1/2- 15) 2 + (0.1/2 -160) 2 - Of! 005 = 5.5 pc, where 
9 denotes the restored beam size, and 1 and c refer to line 
and continuum emission, respectively. 

Fig. [3] presents a close-up of the region around the nu- 
cleus of 3C403, showing the relative distribution of maser 
components I, II, and III. The symbols and error bars 
mark the variance-weighted average locations and uncer- 
tainties of the three maser components derived fitting, 
with the AIPS task JMFIT0, the maser emission in each 
channel with flux density > 5 times the rms noise of 0.6 
mJybeam -1 . Our result indicates an upper limit for the 
offset between different components of 12 ± 4.5 mas. The 
red-shifted emission (components I and II) is seemingly lo- 
cated eastwards w.r.t. the blue-shifted one (component II). 

In addition, we have produced a ~ 3' x 3' map, cover- 
ing the entire radio extent of the galaxy and no additional 
maser spots were detected above the 0.6 mJy level given in 
Sect.H 

4.1.3. CO spectra and millimeter wave continuum 

Motivated by the detection of megamaser emission, we 
searched for the two ground rotational transitions of the 
most common observable molecule, CO. Unlike 22 GHz 
H2O, the CO (1-0) and (2-1) lines predominantly trace 
cool, low-density gas. No clear emission is seen down to 
3(7 noise limits of 13 and 38mJy for the two respective 
transitions (this refers to the unsmoothed spectra with 1 
and 4 MHz channel spacing, see Sect. [2]). We tentatively de- 
tect emission at ~ -300 kms -1 with respect to the sys- 
temic velocity (Sect. [3]). Evans et al. (|2005p also reported 
a non-detection of CO toward 3C403 as part of a survey 
of IRAS radio galaxies. Curiously, also in their spectrum 

5 For a discussion on the accuracy of the parameters derived 
from JMFIT, see e.g. Henkel et al. ((20041 their Sect. 3) 



(their Fig. 1) a highly tentative feature is present at about 
-300 kms -1 . While the offset is opposite to that expected 
from optical emission lines (Sect. [3]) and while evidence for 
such a feature is too weak to seriously question the reces- 
sional velocity adopted throughout this paper, a CO spec- 
trum with higher sensitivity would be desirable. 

Indicated by the baselines, continuum emission was 
also detected. At 109 GHz, the level is approximately 25 
±13mJy (the uncertainty was estimated from the varia- 
tions of the continuum level in individual scans). No con- 
tinuum was detected at 218 GHz. 

4.2. Continuum emission 

4.2.1. 22 GHz VLA emission 

In Fig. 03 we present a naturally weighted VLA A-array 
22 GHz radio continuum image of the nuclear region of 
3C403, produced using maser emission-free channels. The 
source is spatially unresolved with a peak flux density of 
40mJy beam -1 . However, an elongation is apparent in the 
uniformly weighted image of Fig. [6] that might be a weak 
signature of a parsec scale jet visible in the higher resolution 
images presented in Sect. l4.2."2l 

As for the H2O line emission, we have also investigated 
a ~ 3' x 3' area in a search for additional radio continuum 
emission. Apart from the nucleus, among the radio contin- 
uum features (lobes, radio jets, and bright knots) reported 
by BBL and DSL, we have only detected emission at po- 
sition RA J2000 = 19 h 52 m 17.6 s , Dec J2 ooo = 02° 30' 33", 
coincident within the errors with the exceptionally bright 
radio knot F6 (nomenclature according to BBL). The knot 
in our map (Fig. [7]) seems to be resolved into two compo- 
nents aligned almost perpendicular to the large scale radio 
jet. Using a 500 kA taper function, weak emission is also 
detected at the position where two compact structures in 
the eastern hotspot (Fl and F2, following BBL) have been 
previously observed (Fig. [5]). The remaining radio emission 
from 3C403 detected by DSL is either resolved out or is 
below the detection threshold given in Sect.[U 

4.2.2. VLBI maps 

Milliarcsecond resolution 6-cm VLBI maps of the radio nu- 
cleus of 3C403 are shown in Figs. 191 and [11)1 The visibilities 
of 3C 403 demonstrate that the source is resolved at longer 
baselines (Fig.fTTj). The amplitude drops from 70mJy at the 
shortest baselines (10 MA) to 40mJy at 60 MA. 

The uniformly weighted VLBI map of 3C 403 is shown in 
Fig. [9] Notable are extensions to the north-east and south- 
west. We fitted circular Gaussian components to the uv- 
data, and the best fit with three components gave a reduced 
X 2 per degree of freedom of ~ 1.7. Reducing the number of 
components to two or one increases the reduced \ 2 P-d.f to 
2.2 and 2.9, respectively. On the other hand, adding more 
components does not significantly improve the fit. Hence, 
we used a three component model and the results from this 
fit are shown in Table [T] We fitted the three components 
many times with different initial parameters and used the 
scatter in the results to estimate the uncertainties of the 
fit parameters. While the flux densities and angular sizes of 
the three components are not strongly constrained by the 
fit, the position angle is well determined. The orientation of 
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Fig. 5. A naturally weighted 22 GHz VLA radio continuum 
image (resolution: O'.'l, corresponding to ~115pc) of the 
nuclear region of 3C 403, made using channels free of line 
emission. The peak is at 40mJybeam _1 . The contours are 
(-1, 1, 2, 4, 32) x OJSmJybeam" 1 . 
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Fig. 6. A uniformly weighted 22 GHz VLA radio continuum 
image (resolution: 0'.'08, corresponding to ^90 pc) of the 
nuclear region of 3C 403, made using channels free of line 
emission. The peak is at 39mJybeam~ 1 . Contours are (-1, 
1, 2, 4, 32) x 0.75 mJy beam" 1 . 



the pc-scale jet is in very good agreement with the kpc-scale 
jet seen on VLA images. 

Since the telescope in Urumqi (China) could not take 
part in the observations, all information on the long base- 
lines is based on only two telescopes, Hartebeesthoek 
(South Africa) and Shanghai (China). Hence the data on 
long baselines, mainly responsible for the uncertainties in 
the model fit parameters in Table [TJ should be considered 
with some caution. To test the result shown in Fig.[5]we also 
self-calibrated and imaged the data from the short baselines 
only. As before, the amplitude in the visibilities shows a 
slow decrease towards longer baselines. The resulting map 
is shown in Fig.[l0j The extension toward the north-east 
is visible in this map and there might be also a weak jet 



33.8 
33.6 
33.4 



32.8 
32.6 
32.4 
32.2 



O 




O" 



19 52 17.64 17.62 17.60 17.58 17.56 17.54 
RIGHT ASCENSION (J2000) 

Fig. 7. A naturally weighted 22 GHz VLA radio contin- 
uum image (resolution: O'.'l, corresponding to ~115 pc) 
of the bright radio knot (F6, following DSL) in 3C403, 
obtained from channels free of line emission. The peak 
is at 1.5 mJy beam . Contours are (-1, 1, 1.5, 2) x 
0.6mJybeam _1 . 
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Fig. 8. A 500 kA tapered 22 GHz VLA radio continuum 
image (resolution: 0'.'33, corresponding to ~375pc) of the 
compact hotspot structures (Fl and F2, following BBL) in 
3C403, made using channels free of line emission. The peak 
is at 2.5 mJy beam -1 . The contours are (-1, 1, 1.5, 2, 2.5) 
x 0.9 mJy beam -1 . 



component towards the south-west. From this we conclude 
that the weak extensions on mas scales are real and are not 
caused by amplitude calibration errors. 

Assuming the presence of a core and one or two nuclear 
jet components, there are three interpretations: 
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Fig. 9. Full resolution EVN image of the nucleus of 3C 403 
at A=6cm. The contours start at 0.2 mJy and increase 
with a factor of 1.5. The beam size of the observations is 
1.87masx0.94mas at a position angle of 76?9. 
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Fig. 10. Low resolution EVN image of the nucleus of 3C 403 
at A=6cm using only the continental EVN antennas. The 
contours start at 0.35 mJy and increase with a factor of 1.5. 
The beam size is 6.05 masx4.67mas at a position angle of 
67?4. 



Component A is the nucleus. Then we are dealing with 
a one sided jet towards the south-west. 
Component B is the nucleus and components A and C 
are jet and counterjet. 



Radial UV distance along P. A. 65° (10 A) 

Fig. 11. The visibilities of 3C403 clearly show that the 
source is resolved on the longer baselines. 

Table 1. Model fit of the three components A, B, and C to 
the uv-data. S is the flux density, d denotes the separation 
from the phase center. 



Comp. 


S [mjyj 


d [masj 


Size [masj 


P.A. 


A 


9.4±5 


2.41±0.6 


0.69±0.20 


72.5±3 


B 


45.2±10 


0.00 


0.15±0.05 


0.0 


C 


14.7±7 


1.02±0.3 


< 0.04 


-117.1±5 



— Component C is the nucleus and a one sided jet goes 
towards the north-east. 



5. Discussion 

5.1. Evidence for an edge-on torus 
5.1.1. The radio continuum 

To select the most likely of the three scenarios outlined 
in Sect. l4.2."2l and to identify the location of the radio nu- 
cleus, we compare our continuum VLBI maps with the large 
scale structure of 3C 403. The orientation of component A 
(P. A. =73°; Table [T]) is in very good agreement with the 
eastern large scale jet and the position of a bright knot 
(73°) in its lobe (e.g., BBL, DSL). The south-western jet 
(component C), as seen in Fig. [9) seems to bend slightly to 
the west, thus pointing towards the western large scale jet. 
In view of the limited uv-coverage w.r.t. long baselines (see 
Sect. l4.2."2|) . we consider this latter coincidence as tentative. 

The presence of the western and eastern large scale jets 
with similar flux densities suggests that the jets must lie 
close to the plane of the sky. The inner 10 pc of 3C403 
(Figs. l9l and [TO]) show a similar morphology. Thus a scenario 
with B representing the core and A and C marking jet and 
counterjet is by far the most likely one. If we see the jet and 
the counterjet with similar strengths, the orientation of the 
jets must be close to the plane of the sky. This view is sup- 
ported by Kraft et al. (|2005[ hereafter K05). They present 
results from Chandra X-ray observations and describe the 
X-ray spectrum of the nucleus as a superposition of two 
power-law continuum components with a 6.4 keV Fc line. 
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One of these power-law components is heavily absorbing 
(N H ~ 4xl0 23 cm- 2 ). 

If the X-ray photons are predominantly scattered radi- 
ation not passing through the main body of the obscur- 
ing torus, the high column density derived by K05 would 
only be a lower limit. Since no evidence was found for a 
significant amount of dust in observations with the Hubble 
Space Telescope (Martel et al. Q~999) , K05 conclude that the 
absorption is due to material close to the nuclear engine, 
possibly forming a molecular torus. Spectroscopic stud- 
ies at optical wavelengths revealed narrow emission lines 
(Tadhunter et al. 1993), indicating that the broad line re- 
gion is hidden behind an edge-on dusty torus. 

5.1.2. On the absence of Hi and CO 



Morganti et al. (|200ip searched for Hi absorption in 3C 403 
and got an upper limit of < 5xl0 20 cm~ 2 (T sp in = 
100 K). Large differences in Hi and X-ray column densities 
are quite common in Seyfert galaxies (e.g. Gallimore et al. 
1999). The striking discrepancy between X-ray (Sect. 5.1.1) 
and Hi column densities indicates the predominance of ei- 
ther ionized or molecular gas along the line of sight toward 
the nuclear source of 3C 403. While both kinds of environ- 
ment may coexist in a small volume surrounding the central 
engine, the maser emission provides substantial evidence for 
the latter. 

Why are we then not seeing CO? For a quasi-thermally 
excited line, seen in emission, beam dilution may become 
important in an object as distant as 3C403. With the ob- 
served brightness temperature as a function of the intrinsic 
one as well as source and beam solid angle (T£ hs = [T" 1 ' 

X ^torus]/[^torus ^beam 

]) we can determine the small- 
est toroid size, fitorus; we would have detected. With a CO 
line detection threshold of T^ s = 1.8 mK and a 109 GHz 
beam size of 23" (Sect. 2), we obtain fitorus ~ 0'.'04 x 
(500 K/T™') 1 / 2 that corresponds to a linear scale of ~ 50 pc 
x (SOOK/X^ 11 *) 1 / 2 . Elliptical galaxies are characterized by 
high stellar velocity dispersions efficiently heating any cool 
gas and creating an X-ray emitting plasma that is greatly 
confining the volume occupied by molecular clouds (e.g., 
Wikli nd fc Henkel l200Tj ) . The lack of dust emission (Martel 
et al. I1999P supports a small size of the molecular complex 
in the nuclear region and is therefore consistent with our 
non-detection of CO. 

For the absence of deep CO J=1~0 absorption against 
the nuclear radio continuum source (the absorption must 
reach a significant fraction of the continuum level, otherwise 
it would not be detected; see Sect. 14. 1.31 and Fig. Q|, there 
are several possible explanations. (1) There is no molecu- 
lar gas in front of the nuclear source. In view of the col- 
umn density (K05) this is not likely. (2) The torus is partly 
or entirely molecular but has a very small fractional CO 
column density in the J = and 1 states. This may be a 
consequence of a high temperature of the gas, possibly cou- 
pled with an intense radiation field (resulting in excitation 
temperatures of several 100 to ~ 2000 K). (3) The obscur- 
ing clouds may be much smaller than the radio continuum 
source, similar to what was proposed for the much stronger 
radio emitter Cyg A by Barvainis & Antonucci (1994). 
Within this context the size of the H2O maser hotspots 
(Sect. 4.1.1) does not provide any guideline, because 22 GHz 
H2O emission may originate from much higher density gas 



then the ground rotational lines of CO. (4) Alternatively, 
it is also possible that any existing absorption would be 
'swamped' by emission at the same velocity interval. 

5.2. Origin of the H2O maser emission 

5.2.1. Constraints on the application of the standard model 

The main result of our VLA observations is that the two 
22 GHz H2O velocity components of 3C403 do not arise 
from the extended lobes of the system, but are aligned with 
the central radio continuum source within a ~5.5 pc sized 
region. Furthermore, the nuclear torus likely has an edge-on 
orientation (see Sect. 15. ip . Thus an association of the maser 
features with a dense circumnuclear torus/accretion disk is 
a possibility worth discussing. 

To analyze this case, we have to evaluate the statisti- 
cal properties of the sample of known megamaser galaxies. 
Accounting for the sources given in Henkel et al. ([2005), 
Kondratko et al. (|2TjD!)]) . and Zhang et al. (|2"UU5|) . the latter 
containing 64 galaxies with detected H 2 emission beyond 
the Magellanic Clouds, there are a total of 51 known mega- 
maser sources (£h 2 o > 1OL ). Most of these have not yet 
been studied in spatial detail. Omitting 3C403, the sample 
contains two pure jet-masers (apparent interaction between 
a nuclear jet and a dense molecular cloud), while 16 sources 
appear to be disk-masers arising from an accretion disk. 
The latter is either based on detailed observations as in the 
case of NGC 4258 (e.g., Miyoshi et al. [M?5"j) or, a little more 
speculative, on the measured lineshape (see, e.g., NGC 6323 
in Braatz et al. I2004|) . H2O megamaser emission associated 
with a nuclear outflow is only known for one source of the 
sample, the Circinus galaxy. Even here, however, only a 
part of the emission is associated with the outflow, com- 
plementing a nuclear accretion disk viewed approximately 
edge-on (Greenhill et al. I2003p . 

In view of the known sample of megamaser galaxies, an 
association of the maser in 3C403 with a nuclear outflow 
appears to be unlikely. So far all detected 'jet-maser' com- 
ponents are observed at one side w.r.t. the systemic velocity 
(Henkel et al. I2005|) . We should keep in mind, however, that 
3C403 is an (elliptical) FR II and not a (spiral) Seyfert 2 
or LINER galaxy and that "disk-masers" may be identi- 
fied more easily than "jet-masers". Nevertheless, observed 
lincshapes and the relatively large number of disk masers 
superficially hint toward H2O emission from an accretion 
disk as the most plausible interpretation of the H2O lines 
from 3C 403. 

If the emission were part of an accretion disk, at which 
galactocentric radius would it arise? Assuming that the 
two spectral components originate from the tangentially 
seen parts of a nuclear accretion disk seen approximately 
edge-on (i ~ 90°), we obtain (Fig.[T|) a rotation velocity 
of Vrot = Vobs sin -1 i ~100kms _1 . From Bettoni et al. 
([20031 their Table 3) we obtain M B h = 10 8 - 8 M© as the 
mass of the nuclear engine. Combining these two parame- 
ters and assuming the presence of Keplerian rotation (as in 
NGC 4258, see Miyoshi et al. [T9951 Herrnstein et al. IT9991 
2005), this yields with 



R = 0.89 



BH 



M.- 



Vo- 



te sin 



kms 



-1 



D 



Mpc 



an angular distance of 245 mas, corresponding to a galacto- 
centric radius of Rqc ~ 275 pc. This size conflicts with the 
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result of our VLA observations (Sect. E3~l~2")) that, indicat- 
ing that the continuum and the H2O emission are arising 
from the same 5.5 pc (5 mas) sized region, imply a max- 
imal galactocentric radius for the H2O masers of 2.75 pc 
(2.5 mas). 

Is it possible to reconcile this value with the possibil- 
ity that the maser arises in an accretion disk like that in 
NGC4258? A 2.75-pc radius disk rotating at 100 km s" 1 
would imply a black hole mass A/bh = 10 6 8 M Q , a fac- 
tor of 100 lower than that derived by Bettoni et al. 12003 
Alternatively, we might consider a disk inclined by ~6° or 
less. However, this would describe an almost face-on disk, 
a scenario not only different from that of NGC4258 but 
also seemingly contradicting the presence of jets oriented 
approximately along the plane of the sky. 

To summarize, the masers in 3C403 are peculiar in 
three different ways. (1) They are detected at the core of 
a radio galaxy. (2) The two main ultraluminous velocity 
components are highly time variable and (3) the Keplerian 
model seems to fail. If the main velocity components do 
not represent the two tangentially viewed parts of the mas- 
ing disk but instead only one part and the systemic fea- 
tures, this would lead to Vfot~200kms _1 and with A/bh 
= 10 8 - 8 M Q (Bettoni et al. I5005|l to a galactocentric dis- 
tance of 60 mas (70 pc). In this case the disagreement with 
observations is slightly less severe but the linear scale is 
still inconsistent with the results from our VLA A-array 
measurements fSect. l4.lT2"|) . Thus, the accretion disk sce- 
nario requires relatively strict constraints to be applied to 
our case unless we are only viewing the systemic velocity 
features of the putative maser disk (see below) . 

5.2.2. Alternatives to the accretion disk scenario 

Cloud stability considerations: A main difference be- 
tween Seyfert and FR II galaxies (equivalent to radio loud 
type 2 quasars according to the standard unified scheme) 
is the mass of the nuclear engine. While in Seyfert galaxies 
masses range from ~10 6 to a few 10 7 M Q (e.g., Greenhill et 
al. [10001 [TOOT! Herrnstein et al. ITOOOl Henkel et al. 12005)1 . 
masses in radio galaxies and quasars reach 10 9 M Q , with 
3C403 hosting a nucleus that is not far below this latter 
value (Bettoni et al. [217031) . Could clouds in an NGC 4258- 
like disk that is characterized by a particularly large rota- 
tional velocity and small galactocentric radius be stable in 
an environment dominated by such a large central mass? 

The masers in NGC 4258 might arise from well confined 
discrete clumps that are encompassing only a small fraction 
of the entire volume of the warped circumnuclear disk. For a 
roughly spherical clump at distance i?cc from the central 
engine and neglecting magnetic fields (see Modjaz et al. 
121305)) , a density of 

V 21 2nG \R GC J 

is required to reach stability against tidal disruption (Stark 
et al. 11555)1 , G is the gravitational constant. For NGC 4258, 
i.e. i?cc ~ 0.2 pc and V mt ~ 1000 km s" 1 (Herrnstein et 
al. [1999, 2005), the resulting density becomes n(H 2 ) ~ 
5xl0 10 cm -3 . This is close to thermalization (Kylafis et al. 
1991). For commonly accepted collisional excitation this 
either hints at excitation by a two temperature gas (e.g., 
Kylafis et al. 1987) or at emission from turbulent, warm, 



dense molecular debris that does not form self-gravitating 
clouds. 

Most circumnuclear maser disks are characterized by 
smaller rotation velocities and higher galactocentric dis- 
tances, significantly reducing the required density for sta- 
bility against tidal disruption so that here the assump- 
tion of long lived masing clumps giving rise to collision- 
ally excited masers is less of a problem (e.g., n(Ha) <; 
10 8 c m- 3 in NGC 1068; see Greenhill et al. [TITOS) Gallimore 
et al. I200ip . Assuming for a putative accretion disk in 
3C403, where only the systemic spectral features are seen 
(Sect. 5.2.1), R GC = 0.2 pc (as in NGC 4258), the corre- 
sponding minimum density for cloud stability reaches, how- 
ever, 10 12 cm~ 3 , which is prohibitively large for collisionally 
excitated maser emission by a one temperature gas. Could 
then a lack of stable, well ordered clouds explain the rapid 
variability of the maser features in 3C 403? At the maximal 
galactocentric distance permitted by our VLA data, Rqc 
~ 2.75 pc (Sect. 4.1.2), the required density would drop to 
more reasonable JJlO 8 cm~ 3 . 



The jet-maser scenario: While in Seyfert and LINER 
galaxies statistics favor accretion disk over jet-masers (Sect. 
5.1.1), the two main H2O velocity components are (1) quite 
broad and (2) offset by several lOkrns -1 from the nominal 
(NED) systemic velocity of the galaxy. These are the typical 
spectral properties of jet masers (Peck et al.[3t)03 Henkel et 
al. I2005|) . The flares of velocity components I and II (Fig. [2]) 
may arise from shocked regions at the interface between the 
energetic jet material and the molecular gas of the cloud the 
jet is boring through. The time delay between components 
I and II (Fig.[T]) and the apparent red- (component I) and 
blue-shift (component II) with respect to the systemic ve- 
locity may be caused by two clouds at opposite sides of the 
nucleus that have a slightly different galactocentric radius, 
possibly being part of a circumnuclear molecular ring. 



A binary supermassive black hole in 3C403? Aside 
from the mass of the nuclear engine, there may be another 
fundamental difference between 3C 403 and NGC 4258 that 
is worth mentioning. On a large spatial scale, 3C 403 shows 
not the standard two-sided radio jet, but a pair of such 
jets (see Sect. 1). The apparently younger one is oriented 
roughly NE-SW, the more extended, weaker one in the SE- 
NW direction. It is thus one of the rare examples of an 
X-shaped radio galaxy (e.g., Figs. 1 and 2 of Kraft et al. 
2005). This morphology is believed to be caused by the 
coalescence of two supermassive black holes (e.g., Gopal- 
Krishna et al. 120031 Liu I2004p whose interaction modifies 
the inclination of the pre-existing nuclear disk of the dom- 
inant galaxy. Such a scenario may lead to a complexity 
that is much higher than that encountered in the nuclear 
environment of NGC 4258, yielding a puzzling number of 
possible disk- and jet-maser configurations that have to be 
constrained by future measurements. In view of the number 
of free parameters in such an environment, it remains to be 
seen whether the anticorrelation of the main two maser fea- 
tures was a rare accidental event or whether it will provide 
an important clue for a better understanding of the so far 
poorly explored circumnuclear regions of X-shaped radio 
galaxies. 
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6. Conclusions and outlook 

We have investigated the luminous megamaser in the X- 
shaped FR II galaxy 3C 403 using single-dish monitoring of 
the line features and VLA A-array observations. The two 
main maser features show extreme intensity variations dur- 
ing a two-year period. The position of the maser as derived 
from the VLA maps is consistent with a location within a 
few pc of the radio nucleus of the galaxy. 

In addition, we have used the EVN to study the galaxy's 
nuclear radio continuum on parsec scales. A dominant cen- 
tral component and two extensions toward the south-west 
and north-east are identified. Because these parsec-scale 
extensions show position angles that are compatible with 
those of the large scale jets and because these large scale 
jets show similar flux densities, an interpretation in terms 
of a dominant nuclear core and a two-sided jet is plausible. 

No clear sign of CO emission or absorption has been 
detected. If emission and absorption are not accidentally 
canceling each other, this limits the extent of a low density 
molecular gas component to a size of order pc. 

If the mass of the nuclear engine is as large as previously 
estimated and if the nuclear accretion disk is oriented (as 
expected) almost edge-on, an NGC 4258-like accretion disk 
scenario is only viable, if all the detected H2O features are 
systemic, arising from the front or back side of the disk. 
There is, however, no observational evidence for the tan- 
gentially viewed approaching and receding parts of the disk. 
A secular velocity drift of the observed components is also 
not seen. In view of linewidths and deviations from the sys- 
temic velocity, an interpretation in terms of "jet-masers" , 
arising from molecular clouds shocked by the nuclear jet(s), 
is the most plausible one. 

Since 3C 403 is an X-shaped radio galaxy, it may host a 
binary nuclear engine, greatly enlarging the number of po- 
tential scenarios that should be observationally constrained 
by future measurements. These imply single-dish maser 
monitoring over long time spans (1) to study the lifetime of 
the various components, (2) to check whether the anticorre- 
lation of the two main features is a typical or an accidental 
event and (3) to find, possibly, some evidence for a system- 
atic velocity drift. The latter would be a clear sign for the 
presence of a nuclear disk even if the red- and blue-shifted 
parts of this disk remain undetectable (see Wilson et al. 
1995 and Braatz et al. 120031 for such a previously identi- 
fied case). Another very important measurement would be 
the simultaneous observation of the maser features and the 
radio continuum at (sub-)milliarcsecond resolution. This 
could directly show whether the maser features are asso- 
ciated with the nuclear jets or whether we should reject 
this scenario. In any case it would shed light onto the mor- 
phology of the nuclear environment of a prototypical X- 
shaped galaxy. Finally, systematic monitoring of the radio 
continuum, H 2 and X-ray spectra may provide through re- 
verberation mapping important correlations that would be 
essential to assess the physical state and three-dimensional 
morphology of the circumnuclear medium surrounding a 
1O 8 - 8 M AGN. 
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APPENDIX: 4C +02.49 

In order to find suitable phase calibrators for VLBI ob- 
servations, we observed nine bright and compact sources 
from the NRAO VLA Sky Survey (NVSS) with an angu- 
lar separation of <2° with respect to 3C403. The 22 GHz 
(A^1.3cm) observations were carried out with the VLA 
on February 5, 2004. The total observing time was 1.5 
hours yielding ^4 minutes on-source integration time per 
source. The data reduction was performed using standard 
procedures in the Astronomical Image Processing System 
(AIPS). The most suitable source w.r.t. morphology and 
flux density turned out to be 4C+02.49 (see Sect. 2). 

The VLA 1.3 cm observations show that 4C+02.49 is 
a compact double with flux densities of 21.7±0.5 and 
15.3±0.5mJy for the north-eastern and south-western com- 
ponent, respectively (Fig. [12]). A broad band radio spec- 
trum of 4C+02.49 (see the inset in Fig.[T2]) with total inte- 
grated flux densities obtained from the NASA Extragalactic 
Database (NED) and our VLA observation gives a steep 
spectrum with spectral index of -1.07±0.03. Hence, this 
source should be classified as a compact steep spectrum 
(CSS) source (Fanti et al. I2001[) . The tight correlation be- 
tween all data points indicates that there is little variability 
in the total flux. According to the turnover-frequency vs. 
linear size relation in CSS sources (O'Dea & Baum |1997|1 . 
one expects a source size of >10kpc for a turnover fre- 
quency < 100 MHz. This cannot be tested, since, so far, no 
redshift has been derived for this source. 

The two compact components of 4C+02.49 were clearly 
detected with the EVN at A~6cm (Fig.[H]). The flux den- 
sities are 18.4±0.6 (upper left) and 8.8±0.6mJy (lower 
right). Hence, most of the emission that is responsible for 
the total integrated (NED) 6 cm flux of ^220 mJy is re- 
solved out. This is likely the main cause for the seemingly 
inverted spectra when comparing the EVN 6 cm with the 
VLA 1.3 cm flux densities. 
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Fig. 12. VLA map at 1.3 cm (center) and EVN images 
(upper left and lower right) of the two components of 
4C +02.49. Also shown is a broad band radio spectrum (up- 
per right) obtained from NED and our VLA 22 GHz flux. 
The large (blue) crosses indicate flux density measurements 
that had no error bar. 



